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Introduction 


1 . 

In  this  note,  we  add  a  new  wrinkle  to  the  very  old  problem  of 
determining  the  motion  of  a  mass  point  on  a  spring.  We  adopt  a 
general  model  for  the  spring  in  which  the  force  needed  to  compress 
it  to  zero  length  is  infinite.  {Consequently  the  motion  is 
governed  by  a  singular  nonlinear  second-order  ordinary  differen¬ 
tial  equation. )  In  this  setting  we  entertain  the  possibility, 
permitted  by  the  governing  equations,  that  such  a  total 
compression  is  actually  attained.  This  total  compression  corres¬ 
ponds  to  a  kind  of  shock.  We  then  extract  from  the  governing 
equations  all  the  illumination  they  can  shed  on  the  physical 
behavior.  Our  problem,  which  presents  novel  features  for  ordinary 
differential  equations,  captures  in  microcosm  deep  and  unresolved 
issues  involving  shocks  and  their  suppression,  which  arise  in  the 
study  of  quasilinear  hyperbolic  and  parabolic  partial  differential 
equations.  We  comment  briefly  on  these  issues  in  Section  8. 

Notation.  t  denotes  the  time.  Ordinary  derivatives  with  respect 
to  it  are  denoted  by  superposed  dots.  All  other  ordinary 
derivatives  are  denoted  by  primes.  We  denote  partial  derivatives 
by  subscripts.  Certain  symbols  such  x  appear  in  two  roles:  as 
a  real-valued  argument  of  a  function  and  in  an  expression  x(t) 
for  the  value  of  a  real-valued  function  of  t.  For  the  sake  of 
precision  we  denote  the  function  with  values  x(t)  by  x(*)* 

2.  Formulation  of  the  Governing  Equations 

A  spring  (e.g.,  in  the  form  of  a  steel  helix  or  a  rubber 
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cylinder)  is  confined  to  a  horizontal  groove  along  the  x-axis. 

The  "left"  end  of  the  spring  is  fixed  to  a  wall  at  x  =  -  1.  A 
unit  mass  point  is  attached  to  the  "right"  end.  Without  loss  of 
generality  we  assume  that  the  natural  (unstretched)  length  of  the 
spring  is  1 .  We  let  x(t)  denote  the  position  of  the  mass  point. 
(x(t)  is  also  the  displacement  of  the  mass  point  from  its 
equilibrium  state.)  See  Fig.  1. 

We  assume  that  the  force  exerted  at  time  t  by  the  spring  on 
the  mass  point  depends  only  on  the  position  x(t)  and  the  velocity 
x(t)  of  the  mass  point.  We  denote  this  force  by  -n(x( t ) ,x( t ) ) . 
The  function  n  is  assumed  to  be  given.  Let  us  set 

(1)  n(x,0)  =  ?>'(x)  with  <p'{0)  =  0,  <p{0)  =  0. 

Then  n  has  the  form 

(2)  n(x,y)  =  <p'  (x)  +  yg(x,y)  with  g(x,0)  =  0. 

Here  v  is  a  non-negative  parameter.  <p'(x)  is  the  tensile 
force  exerted  on  the  mass  point  when  it  has  position  x  and  is 
not  moving.  It  effectively  describes  the  elastic  properties  of 
the  spring.  g  accounts  for  forces  that  come  into  play  when  the 
mass  point  moves.  These  are  the  internal  frictional  or  viscous 
forces  of  the  spring.  When  the  mass  point  moves  both  the  groove 
and  the  ambient  air  can  exert  frictional  forces  on  the  mass  point. 
Such  forces  can  be  absorbed  in  g.  For  our  purposes,  however,  it 
is  conceptually  advantageous  to  regard  these  external  forces  as 
absent  (in  consequence  of  a  perfect  lubrication  of  the  groove  and 
of  the  placement  of  the  system  in  a  vacuum) .  We  assume  that  no 
other  forces  act  on  the  spring.  Then  the  equation  of  motion  for 
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the  mass  point  is 


(3)  x(t)  +  <p'(x(  t))  +  pg{x( t) ,x(t) )  =  0. 

For  simplicity  we  assume  that  n  is  continuously 
differentiable.  We  require  that 

(4)  n  (x,y)  >  0  Vx>-1  and  Vy 

so  that  an  increase  in  the  length  of  the  spring  is  accompanied  by 
a  corresponding  increase  in  the  tensile  force.  We  adopt  growth 
conditions  compatible  with  (4):  For  each  fixed  y 

(5)  n(x,y )  — >  as  x  >  {jfj- 

so  that  an  infinite  tensile  force  is  needed  to  stretch  the  spring 
to  infinite  length  and  that  a  negatively  infinite  tensile  force 
(i.e.,  an  infinite  compressive  force)  is  needed  to  compress  the 
spring  to  zero  length.  Conditions  (4)  and  (5)  yield  corresponding 
conditions  on  <p' .  We  illustrate  <p  in  Fig.  2.  (<p  is  called  the 
stored  energy  of  the  spring).  That  <p(x)  need  not  approach  ® 
as  x  — ►  -  1  is  crucial  in  our  development. 

We  assume  that 

(6)  ny(x,y)  >  0. 

Thus  an  increase  in  the  velocity  of  the  end  of  the  spring  is 
accompanied  by  a  corresponding  increase  in  the  tensile  force.  As 
we  shall  see,  this  requirement  ensures  that  the  effect  of  g  is 
to  dissipate  energy.  When  v  *  0,  we  accordingly  say  that  the 
motion- is  damped .  We  adopt  a  strong  version  of  (6)  that  is 
required  to  hold  when  the  spring  is  under  compression  (x  <  0)  and 
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is  being  compressed  (x  <  0):  There  is  a  positive-valued  function 
¥>'  such  that 

(7)  g(x,y)  s  ?'(x)y  for  -1  <  x  s  0,  y  s  0. 

(Note  that  if  (7)  is  an  equality,  then  (4)  holds  for  all  x  in 
(-1,0]  and  for  all  y  s  0  if  and  only  if  y>"  (x)  so.)  We  take 
y>(  0)  =  0.  The  viscous  stress  g  is  allowed  to  depend  upon  x 
because  the  internal  frictional  force  when  the  spring  is  highly 
compressed  may  well  differ  from  that  when  the  spring  is  under 
tension. 

3.  Analysis  of  the  Differential  Equation  for  the  Undamped  Motion 

We  begin  our  study  of  (3)  by  first  studying  it  for  v  =  0. 

In  this  case  we  multiply  (3)  by  x(t)  and  integrate  the  resulting 
equation  to  get  the  energy  equation: 

(8)  §x(t)2  +  *(x(t)}  =  |k(0)2  +  *>(x(0)). 

This  integral  is  assumed  to  hold  only  on  time  intervals  on  which 
x(‘)  is  twice  differentiable  and  satisfies  x(t)  >  -  1  for  all 
t's  in  the  interval. 

We  often  set  y(t)  =  x(t).  The  phase  portrait  of  (3) 
consists  of  the  family  of  curves 

(9)  5(E):  |  y2  +  *(x)  =  E(const) 

in  the  (x,y)- plane.  Equation  (8)  says  that  a  phase-plane  trajectory 
(X(*)#x(*))  Of  the  equation  is  confined  to  a  single  curve  of  the 
form  (9)  on  any  time  interval  on  which  x(«)  is  twice  differen¬ 
tiable.  (Since  <*>(x)  — >  <»  as  x  — ►  oo,  equation  (8)  implies  that 
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there  is  a  number  C  (depending  on  initial  condition)  such  that 
|x(t)|  <  C,  x(t)  <  C.  Standard  continuation  results  for  ordinary 
differential  equations  (cf.  [2,  Chap.l.])  then  ensure  that  the 
solution  of  (3)  exists  as  long  as  x(t)  >  -  1.) 

If  $>{x)  — >  oo  as  x  — »  -  1,  the  phase  portrait  is  given  by 
Fig.  3,  in  which  no  trajectory  flirts  with  the  wall  x  =  -  1 
forming  the  edge  of  its  universe.  But  if  <p  has  a  finite  limit 
as  x  — >  -  1,  then  the  phase  portrait  has  an  entirely  different 
character  illustrated  in  Fig. 4:  For  large  enough  E  the 
trajectories  actually  touch  the  wall  tangentially.  If  this 
happens  to  a  solution,  i.e.,  if  there  is  a  r  such  that 
x(r)  =  -  1,  we  say  that  the  spring  suffers  a  shock.  It  is  not  at 
all  clear  what  should  happen  to  x(’)  after  the  first  time  t 
that  x(t)  =  -  1.  Our  goal  is  to  deduce  the  fate  of  x(-)  i.e., 
to  deduce  the  shock  structure,  on  the  basis  of  sound  mathematical 
arguments  and  physical  hypotheses. 

Let  us  note  that  the  intersection  of  the  region  enclosed  by 
(9)  and  the  half-plane  {(x,y):  x  >  -  1)  is  convex  because  <p  is 
a  convex  function.  Since  <p  has  its  minimum  at  x  =  0,  the  curve 
(9)  has  its  greatest  separation  in  the  y-direction  at  x  =  0. 

4.  The  Equation  for  Damped  Motion 

We  now  examine  the  equation  (3)  (with  v  >  0)  for  damped 
motion,  which  is  equivalent  to  the  system 

(10a, b)  •  x(t)  =  y(t),  y ( t )  =  -  ?'(x(t))  -  vg(x( t) ,y(t) ) 

of  first-order  equations.  Equation  (10a)  says  that  trajectories 
of  the  phase  portrait  move  right  when  y  >  0  and  move  left  when 
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y  <  0.  The  horizontal  isocline  of  (10)  is  the  locus  of  points  at 
which  trajectories  are  horizontal,  i.e.,  the  set  of  (x,y)  such 
that 

(11)  <p'(x)  +  yg(x,y)  =  0. 

Conditions  (4),  (5)  support  a  global  implicit  function  theorem 
(based  on  the  intermediate  value  theorem)  ensuring  that  (11)  has  a 
unique  solution: 

(12)  x  =  h(y;p) 

where  h  is  a  continuously  differentiable  function  with 

(13)  h(0;i>)  =  0,  hy(y;t>)  s  0,  h(y;l>)  >  -  1,  h(y;0)  =  0. 

The  vertical  isocline  of  (10)  is  the  locus  of  points  at  which 
trajectories  are  vertical;  it  is  the  x-axis  as  a  consequence  of 
(10a) . 

Let  us  multiply  (3)  by  x(t)  obtaining 

(14)  x(t)2  +  <p(x(t))j  =  -  yg(x( t)  ,x(t)  )x( t)  . 

For  v  >  0,  the  properties  of  g  ensure  that  the  right-hand  side 
of  (14)  is  negative  when  x(t)  *  0,  i.e.,  when  a  phase-plane 
trajectory  is  not  crossing  the  x-axis  (cf.  (10a)).  It  follows 
from  (14)  that  trajectories  pierce  inward  through  the  level  curves 
(9).  (This  fact  supports  an  existence  theorem  like  that  discussed 
in  the  remarks  following  (9).)  Equation  (14)  also  says  that  the 
energy,  which  is  the  term  in  braces  on  the  left  side  of  (14),  is 
decreasing  for  v  >  0.  It  is  this  fact  that  justifies 
characterizing  g  as  dissipative. 
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5. 


Shock  Absorbers 


The  dissepativity  embodied  by  (6)  ensures  that  any  trajectory 
(10)  originating  within  or  on  the  distinguished  curve  £(<*>(- 1)) 
of  (9),  which  is  tangent  to  the  line  x  =  -  1,  must  ultimately 
wind  down  to  the  origin.  (This  fact  is  a  standard  result  from  the 
theory  of  ordinary  differential  equations.)  But  it  is  not  clear 
whether  a  typical  trajectory  starting  outside  6(<p(- 1))  will 
intersect  the  line  x  =  -  1,  thereby  generating  a  shock,  or  will 
safely  negotiate  this  dangerous  edge  and  ultimately  bring  the 
system  to  rest  at  (x,y)  =  (0,0). 

Suppose  that  <p  has  a  finite  limit  as  x  — >  -  1  (or 
equivalently,  that  <p'  is  integrable  on  (0,1]).  Then  g  is 
called  a  shock  absorber  if  for  each  v  >  0  there  is  no  trajectory 
starting  to  the  right  of  the  line  x  =  -  1  that  reaches  this 
line.  We  determine  conditions  that  ensure  that  g  is  a  shock 
absorber  and  then  show  how  shock  absorbers  can  be  used  to  develop 
a  shock  structure  for  the  undamped  equation. 

We  integrate  (3)  from  0  to  t  obtaining 

ft  •  r* 

(15)  0  =  x(t)-x(0)  +  v  g(x(r )  ,x(r  )dr  +  <p'(x(r))dr. 

J  0  J  0 

We  study  (3)  and  (15)  on  a  time  interval  [0,t]  .  for  which 
(x(-),x(*))  lies  entirely  in  the  quadrant  Q  s  { (x,y) :-l<x<0 ,y<0) . 
Without  loss  of  generality  we  take  initial  conditions  with 
x(0)  <  0,  x(0)  <  0.  Every  trajectory  in  Q  may  be  regarded  as 
originating  from  such  an  initial  point.  Condition  (7)  implies 
that 
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r1  •  r 

(16)  0  <  X  ( t )  -x(  0 )  +  v  v>'(x(T)  )x(r)dT  +  <p'(x(  T))dr 

J  a  ^  n 


■t  .  r* 

p'(x(r)  )x(r)dT  + 

0  0 

for  trajectories  in  Q.  Since  <p'(0)  =  0,  condition  (4)  implies 

that  <p'(x)  <  0  for  x  <  0.  Thus  inequality  (16)  immediately 

yields 

(17)  vtp(x(  t))  >  x  (  0 )  -x  ( t )  +  v  ¥><x(0))  -  *>'(x(r))dT 

J  0 

>  x( 0 )  +  vv>(x(0) ) 


for  trajectories  in  Q. 

Now  x(0)  <  0  in  Q.  If  y>(x)  — -  co  as  x  *  -  1 ,  then 


(18) 


x  ( t )  *  0))  +  x(0  )/u)  >  -  1 


for  v  >  0 


so  that  g  is  a  shock  absorber.  If  y>(x)  has  a  finite  limit  a 
as  x  “  -  1,  then  there  are  initial  data,  including  those  with 
x(0)  <  vo, ,  for  which  shocks  must  occur.  See  Fig.  5.  (Note  what 
happens  as  v  “  0 . ) 

The  physical  significance  of  the  requirement  that  rp  (x)  - 
as  x  *  -  1  is  that  internal  friction  increases  without  bound  as 
the  spring  is  subjected  to  increasingly  severe  compression. 


-  00 


6.  Shock  Structure 

Suppose  that  the  elastic  energy  <p(x)  has  a  finite  limit  as 

x  _ »  -  i,  in  which  case  the  undamped  problem  has  phase  portrait 

Fig. 4.  Let  g  be  a  shock  absorber.  We  adopt  the  view  that  the 
physically  natural  behavior  for  solutions  of  the  undamped  problem 
in  the ’presence  of  shocks  is  that  dictated  by  the  nonuniform  limit 
as  v  *  0  for  the  problem  damped  with  a  shock  absorber. 


9 


To  determine  this  behavior,  let  us  assume  that  v  is  small. 

Then  standard  perturbation  results  of  ordinary  differential 
equations  (cf.  [2,  Chap. 2],  e.g.)  tell  us  that  away  from  the  wall 
x  =  -  1  (where  the  equation  ceases  to  be  regular) ,  the  ) 

trajectories  of  (3)  are  close  to  those  of  the  undamped  equation. 

(I.e.,  on  any  compact  subset  of  {(x,y):  -1  <  x}  the  trajectories 
of  (30)  uniformly  approach  those  of  the  undamped  equations  as 
v  »  0.  ) 

We  illustrate  a  typical  trajectory  in  Fig.  6  with  large 
initial  data.  It  stays  close  to  the  level  curve  <?(E)  until  it 
approaches  the  wall.  Since  g  is  a  shock  absorber,  the 
trajectory  cannot  touch  the  wall.  Its  motion  in  Q  satisfies 
x(*)  <  0,  y(*)  >0  so  that  it  must  move  nearly  parallel  to  the 
wall  in  Q.  It  is  therefore  constrained  to  enter  the  second 
quadrant.  But  this  trajectory  must  pierce  inwardly  through  every 
level  curve  (9).  Therefore,  in  the  second  quadrant  it  cannot 
cross  the  level  curve  6(<p(~  1)),  which  is  tangent  to  x  =  -  1. 

It  is  therefore  subsequently  confined  to  remain  within  this  level 
curve,  slowly  winding  down  to  the  equilibrium  state  (x,y)  =  (0,0). 

The  behavior  of  the  phase  portrait  in  the  limit  as  v  “  0  is 
nearly  obvious.  The  limiting  trajectory  corresponding  to  that  of 
Fig.  6  is  given  in  Fig.  7.  It  consecutively  passes  through 
points  0,  1  ,  1  ,  2,  3,  4  with  the  cycle  1  ,  2,  3,  4  repeated 
periodically.  The  only  question  facing  us  to  determine  what  is 
happening  on  the  segment  1  ,1+.  To  answer  this  we  revert  to 
the  study  of  (3)  when  v  is  small.  Let  e  be  a  small  positive 
number.  We  study  the  phase  portrait  of  (3)  or  (10)  in  the  strip 
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{(x,y):  -l<x<-l+£,y<0).  Here  y(t)  >  0.  Thus  let 
y(*),  defined  on  (t^tg),  correspond  to  a  trajectory  lying  in 
this  strip.  Let  t ( • )  denote  its  inverse,  defined  on  (y^yg), 
where  y1  =  =  V(t2)-  The  time  lapse  in  traversing  this 

portion  of  the  trajectory  is 

y2  .  Y2  Y2 

(19)  t  -t  =  f  t'(y)dy=f  ^-5^ —  =-j  - = - ^ - - - , 

JYX  Y1y(t(y))  y^'  (x(  t{y)  )  )+vg(x(t(y)  )  ,y) 


the  last  equality  coming  from  (10b).  Since  x 
from  (19)  that 


(20) 


Y2"Y1 

t2  tl  S  *>'(-!+£  )  ' 


< 


-1 


+  £  , 


we  find 


As  £  *  0,  t2-t i  »  0.  We  have  seen  in  Fig.  6  that  as  v  “  0 
the  trajectory  must  get  closer  and  closer  to  x  =  -  1  in  quadrant 
Q.  In  virtue  of  (20)  we  therefore  conclude  that  the  appropriate 
time  lapse  for  the  trajectory  to  go  from  1  to  1+  in  Figure  8 
is  the  limit  of  the  time  lapse  needed  to  traverse  the  nearly 
vertical  part  of  Fig.  6;  this  limit  is  0.  Thus  we  may 
characterize  the  shock  structure  of  the  undamped  problem  as  being 
embodied  in  the  phase  portrait  of  Fig.  7  with  the  jump  from  l- 
to  1  occurring  instantaneously.  Thus  x(*)  is  discontinuous. 
(A  discontinuity  in  the  velocity  field  is  a  difining  property  of 
shock  in  the  theory  of  partial  differential  equations  of 
mechanics.  Our  notion  of  shock  is  consequently  consistent  with 
this  usage.)  We  now  plot  in  Figs.  8  the  position  x(*)  and 
velocity  x(*)  in  light  of  this  interpretation  of  Fig.  7. 
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This  development  has  some  interesting  consequences.  Suppose 
that  <p  is  integrable  on  (-1,0)  and  that  g  is  a  shock 
absorber.  Let  us  integrate  (3)  with  respect  to  t  over  (t^tg) 
obtaining  the  Impulse-Momentum  Principle 


which  is  to  hold  for  (almost)  all  intervals  (t^tg).  The  left 
side  of  (21)  is  the  difference  in  the  linear  momentum  (since  the 
mass  is  unity)  and  the  right  side,  the  time  integral  of  the  force, 
is  the  impulse .  The  Impulse  -  Momentum  Principle  generalizes 
Newton's  Second  Law  of  Motion  to  handle  cases  in  which  the 
velocity  x(*)  is  not  a  differentiable  function  of  t.  More 
generally,  we  could  regard  the  impulse  over  (t^tg),  not  as  an 
integral  of  a  prescribed  force  function,  but  rather  as  a 
prescribed  functional  of  x(*)* 

Suppose  that  the  situation  illustrated  in  Fig. 6  prevails. 

Let  v  “  0.  Then  the  trajectory  of  Figure  7  approaches  the  line 
x  =  -  1  and  in  the  limit  touches  it  along  1  ,1+  of  Fig. 7  for 
an  instant  r  of  time.  Let  (t^tg)  3  T  and  ^et  tg“T. 

Since  x  suffers  a  jump  at  r  in  Fig. 7,  the  formal  limit  of 
(21)  as  v  “  0  cannot  hold.  What  happens  is  that  the  integrand 
on  the  right  hand  side  of  (21)  converges  in  the  sense  of 
distributions  to  a  Dirac  delta.  This  means  that  it  would  be  a 
fruitless  exercise  to  seek  a  solution  of  (21)  with  v  =  0  in  the 
class  of  functions  x(*)  for  which  *>'(x(*))  is  integrable.  In 
this  limiting  case  the  right-hand  side  of  (21)  should  be 
interpreted  as  the  value  of  a  functional  of  x( • ) ,  depending  on 
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Let  us  also  note  that  the  energy  drops  across  r.  Indeed, 
the  impulse  functional  can  be  expressed  in  terms  of  the  energies 
at  r-  and  r+.  (cf- (8) ) -  Thus  in  the  limit  as  v  approaches 
0,  the  effects  of  dissipation  become  concentrated  at  the  shock. 


7.  Example.  Quadratic  Damping. 

We  can  illuminate  the  nature  of  the  inequalities  (6),  (7)  by 
examining  a  situation  in  which  they  are  not  satisfied,  that  of 
quadratic  damping: 

(22)  g(x,y)  =  |  V'(x)y[y|,  y>'(x)  2  0. 


On  any  t-interval  for  which  x(t)  <  0,  we  can  let  t(-)  be 

the  inverse  of  x{  ■  )  and  set  v(x)  =  x(  t  (x)  )  .  Then 
x(t(x))=*v(x)v'(x)  and  (3)  reduces  to 


(23) 


dx  L2 


<xr  +  *>'(x)  - 


^'(x)v(x)2  = 


0  where  x(t)  <  0 


(The  sign  of  the  last  term  on  the  left-hand  side  of  (23)  changes 
•  1  2 

where  x(t)  >  0.)  Let  w  =  ^  v  .  Since  (23)  is  a  linear 
equation  for  w(*),  its  solution  is 


(24)  w(x)  =  d? 


We  assume  that  <p(x)  has  a  finite  limit  as  x  — >  -  1.  Noting 
that  the  independent  variable  In  (24)  is  x,  we  refer  to  Figure  6 
to  deduce  that  (22)  is  a  shock  absorber  if  for  each  w(0)  >  0  and 
for  each  v  in  (0,1]  there  is  an  x  in  (-1,0)  for  which 
w(x)  =  0.  In  this  case,  (24)  implies  that  there  is  an  x  in 
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(-1,0)  such  that 


(25)  w(0)  =  f  <p'(Z  )  e~Uy>iK  }dZ  . 

J  0 

If  y(x)  has  a  finite  limit  as  x  — »  -  1,  then  the  integrability 
of  *'  ensures  that  the  right-hand  side  of  (25)  has  an  upper 
bound  independent  of  i>  for  v  s  1.  Thus  (22)  cannot  be  a  shock 
absorber.  If  y>(x)  — >-co  as  x  — >  -  1,  the  situation  is  far 
more  delicate  than  that  for  (6),  (7),  discussed  in  Section  5. 

E.g.,  if  *(x)  =  -Cn(l+x),  then  the  right-hand  side  of  (25) 
reduces  to 

rx 

(26)  <t>'{Z)(  i+O  az. 

J  o 

*'(x)  >  -  C(l+U”^, 
v  <  1  -  p.  Thus 
(22)  fails  to  be 
-  (l+x)-0  with 
right-hand  side  of 


(l+x)"01 

>  [-<*>'(-1/2)  ]a_1  J  tl+^?7  +  *  '  •  +  (^r?)k/k!  ]77_(a+1)/ad7? 

(1/2)"® 


If  there  are  numbers  C>0,  p  in  (0,1)  such  that 
then  (26)  is  bounded  above  for  -1  <  x  <  0  when 
even  though  y>(x)  — >-<»  as  x  — *  -  1 ,  equation 
a  shock  absorber.  On  the  other  hand,  if  y>  (x)  = 
a  >  0  for  x  <  -  1/2,  then  for  x  <  -  1/2  the 
(25)  exceeds 

-1/2 

-a 


[-*'(-1/2)  ]  |  e‘ 


v{l-Z  ) 


d  Z 


x 


(l+x) 


(27)  =  [-*'(-1/2) 3a"1  J 


(  1/2) 


-a 

e 

-a 


W  ^-(a+D/a 


dr) 
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1 ,  the 


( i+x) ~a 

>  yk[-<*>'(-l/2)](ak!)_1  j  !7k_(a+1)/c'dr? 

(l/2)_Ct 

where  k  is  an  integer  exceeding  (a+l)/a.  Then  as  x  »  - 
rightmost  term  of  (27)  approaches  ao  for  each  v  >  0.  Thus  (22) 
with  this  class  of  rp's  are  shock  absorbers. 

8.  Conclusion 

Our  treatment  in  Sections  4-6  produced  a  doctrine  for  con¬ 
structing  solutions  of  the  undamped  equation  near  a  time  r  at  which 
x(t)  =  -  l.  One  interpretation  of  this  doctrine  is  that  it  yields 
a  criterion  for  selecting  a  unique  way  to  continue  the  solution 
after  the  shock.  Without  such  a  well-founded  doctrine  we  could 
have  invoked  a  spurious  argument  based  on  the  conservation  of  energy 
to  infer  that  a  trajectory  arriving  at  (x,x)  =  (-l,-y0)  should 
leave  along  the  tail  of  the  same  trajectory  at  (x,x)  =  (-l,yQ). 

One  can  contemplate  problems,  analogous  to  the  one  we 
treated,  for  the  longitudinal  motion  of  a  viscoelastic  rod.  This 
motion  is  governed  by  a  partial  differential  equation  of  the  form 

(28)  p(s)utt(s,t)  =  n(us(s,t) ,ust(s,t) )s 

=  [<p'(us(s,t)  )  +  i>g(ug(s,t)  ,ust(s,  t)  )  ]g 

where  p  is  the  mass  density  of  the  rod  per  unit  reference 
length,  u(s,t)  is  the  displacement  of  the  material  point  s  at 
time  t,.  and  n  has  the  form  (2).  Pioneering  analyses  of  the 
existence  of  solutions  for  initial-boundary  value  problems  for  (28) 
were  carried  out  by  Greenberg,  MacCamy,  &  Mizel  [6]  for 
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(29) 


n(x,y)  =  <p'(x)  +  fy,  v  >  0, 


by  Kanel'  [8]  (and  later  independently  by  MacCamy  [19])  when 
there  is  a  number  C  >  0  such  that 

(30)  n(x,y)  =  <p'(x)  +  ?'(x)y,  y>‘  (x)  ;>  C-1, 

and  by  Dafermos  [3]  when  there  is  a  number  C  >  0  such  that 

(31)  ny  ^  C-1  ,  (nx)2  s  Cny. 

It  is  clear  that  if  <p‘  has  the  form  show  in  Fig.  2,  then 
(29)  satisfies  ( 4 ) — ( 7 ) .  Obviously  (30)  satisfies  (6), (7).  But  if 
V'(x)  is  not  constant,  then  (4)  cannot  hold  for  all  y.  Under 
mild  additional  conditions  on  n  it  can  be  shown  [1]  that  the 
second  condition  of  (31)  is  incompatible  with  the  second  limit  of 
(5)  . 

These  observations  do  not  deprive  (30)  and  (31)  of  their 

physical  significance.  Equation  (28)  subject  to  (30)  is  the 

material  (or  Lagrangian)  form  of  the  one-dimensional  Navier-Stokes 

equations  for  the  flow  of  a  compressible  viscous  fluid,  provided 

that  the  first  limit  of  (5)  is  replaced  with  the  inequality 

n(x,y)  <  0.  Here  u  +1  is  the  specific  volume  and  -n  is  the 

s 

pressure.  Equation  (28)  subject  to  (31)  furnishes  a  perfectly 
good  model  for  one-dimensional  shearing  motion  of  a  viscoelastic 
body . 

To  the  best  of  my  knowledge,  there  is  no  global  existence 
theory,  for  initial-boundary  value  problems  for  (28)  for  a  general 
class  of  functions  n  satisfying  (4)-(6).  The  difficulty  lies  in 
the  treatment  of  the  possibility  that  u(s,t)  =  -  1  for  some  s,t 
and  thus  reflects  a  very  challenging  version  of  the  analogous 
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obstacle  overcome  in  Sections  4-6.  (In  technical  terms,  the 
uniform  parabolicity  of  (28)  inherent  in  the  assumption  that 
n^  £  C-1,  is  insufficient  to  overcome  the  extreme  hyperbolicity 
inherent  in  the  second  limit  of  (5).}  I  conjecture  that  an 
appropriate  existence  theory  can  be  fashioned  on  the  basis  of  (7). 

The  limiting  process  of  Section  6  represents  one  particularly 
attractive  way  to  develop  a  doctrine  for  handling  shocks.  Hopf 
[7]  was  able  to  effect  a  very  elegant  treatment  of  the  limiting 
process  for  a  special  partial  differential  equation,  but  the  few 
extensions  of  it  to  more  general  equations,  by  DiPerna  [4,5],  have 
required  formidable  exercises  in  analysis.  In  all  these  cases  the 
dissipative  terms  have  the  special  character  of  (30). 
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Captl 

Fig .  1 
Fig. 2 


Fig. 3 

Fig. 4 

Fig. 5 


Fig .  6 


Geometry  of  the  spring-mass  system. 

Typical  stored  energy  functions  <p »  The  actual  graphs 
shows  are  for 

9>1(x)  =  |(l+x)3  +  (1+x)-1  -  |  , 

V>  2  ( x )  =  (l+x)3/2  -  3(l+x)1/2  +  2. 

Typical  phase  portrait  for  the  undamped  equation  when 
<p(x)  — >oo  as  x  “  -  1 .  That  plotted  here  is  for  <p 
introduced  in  the  caption  of  Fig. 2. 

Typical  phase  portrait  for  the  undamped  equation  when 
<P(-1)  <oo.  That  plotted  here  is  for  <p  introduced  in 
the  caption  of  Fig. 2. 

A  function  ip  that  generates  a  shock  absorber.  That 
plotted  here  has  the  formula 

ip{x)  =  (1+x)  -  (1+x)  1 . 

The  dashed  line  corresponds  to  a  ip  with  >  -  oo. 

Typical  phase  portrait  of  the  damped  system  10  when  g 
is  a  shock  absorber.  The  trajectory  shown  is  computed 
for  <p2  of  Fig.  2,  for  v  =  0.1  and  for 

ip' [x)  y  for  y  <  0 
y  for  y  2  0 

where  ip  has  the  form  given  in  the  caption  of  Fig.  5. 
(That  g  is  not  continuously  differentiable  has  no 
appreciable  effect  on  the  phase  portrait.  This  g  can 


g(x,y) 


■{ 
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Fig. 7 
Fig. 8 


be  readily  approximated  by  a  function  with  any 
prescribed  number  of  derivations  for  which  (4)  -  (7) 
hold.  ) 

The  limit  as  v  *  0  of  the  trajectory  of  Fig. 6. 

The  displacement  x(  •  )  and  the  discontinuous  velocity 
v(*)  for  the  trajectory  of  Fig. 7. 
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